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Matsumoto and Yoshimura have recently argued that the 
number density of heavy particles in a thermal bath is not 
necessarily Boltzmann-suppressed for T -C M, as power law 
corrections may emerge at higher orders in perturbation the- 
ory. This fact might have important implications on the de- 
termination of WIMP relic densities. On the other hand, the 
definition of number densities in a interacting theory is not 
a straightforward procedure. It usually requires renormal- 
ization of composite operators and operator mixing, which 
obscure the physical interpretation of the computed thermal 
average. We propose a new definition for the thermal aver- 
age of a composite operator, which does not require any new 
renormalization counterterm and is thus free from such ambi- 
guities. Applying this definition to the annihilation model of 
Matsumoto and Yoshimura we find that it gives number den- 
sities which are Boltzmann-suppressed at any order in per- 
turbation theory. We discuss also heavy particles which are 
unstable already at T = 0, showing that power law corrections 
do in general emerge in this case. 
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I. INTRODUCTION 

The computation of number densities of cosmological relics, 
such as weakly interacting massive particles (WIMPs), neu- 
trinos, baryon and lepton number, is usually accomplished by 
using a thermally averaged Boltmzann equation |lU2j. The 
picture emerging in this context is that of a sudden freeze- 
out: the particle number density follows the equilibrium value 
N(T) for temperatures higher than the freeze-out value, Tf, 
below which the annihilation is frozen. In the WIMP sce- 
nario for dark matter, the freeze-out temperature turns out 
to be typically much lower than the WIMP mass, then one 
expects a Boltzmann-suppression factor exp(— M/Tt) in the 
final WIMP number density. In a series of papers |H], Mat- 
sumoto and Yoshimura (MY) have challenged the above con- 
ventional conclusion, claiming that two-loop corrections to 
the number density exhibit only power-law suppression in 
Tf/M, thus dominating over the Boltzmann-suppressed tree- 
level contribution if Tf is low enough. If confirmed, such a 
finding would imply that the present constraints on dark mat- 
ter models are actually underestimated and should be care- 
fully reconsidered. To be definite, MY introduced a toy model 
of two real scalar fields with lagrangian 
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with M ^> m. The following relations were assumed among 
the coupling constants 



IAJ < A^ 



|A| < IAJ < 1, 



so that the light % particles act as an efficient heat bath for 
the heavy <j>'s. They then considered the quadratic part of 
the Hamiltonian for <h, 



^ = i0 2 + i(V<« 2 + iM 2 



(2) 



and, for temperatures T such that m <T <C M they defined 
the number density of <j> particles as 



N,u 



{H^ + counterterms)T 
M 



(3) 



where MY define the thermal average for an operator A as 
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tre- H / T 
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H being the total Hamiltonian. At tree-level, the well known 
Boltzmann-suppressed contribution is obtained 
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N w = l^-\ e~ M/T . 



The power-law contribution arises at two-loops and is given 
by 



Af > = cA : 
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which, as anticipated, dominates over the tree-level contribu- 
tion at low temperatures. The interpretation of (0) as the 
physical number density of heavy <f) particles has been later 
questioned by Singh and Srednicki (4J5|. First, they pointed 
out that the splitting of the full Hamiltonian in two parts - 
one for the "system" of <j> particles, and the other for the "en- 
vironment" of x's - is problematic, as the interaction energy 
is not properly taken into account. It was actually shown 
to be of the same order of the "system" energy in an ex- 
actly solvable Caldeira-Leggett model. Second, Srednicki |j] 
considered a modification of the model in eq. (|l|) by taking 
complex (j> fields charged under a (7(1) symmetry. Then, he 
showed that from the fluctuations of the conserved charge it 
is possible to define a number density which exhibits the ex- 
pected Boltzmann-suppression at low temperature. Extrap- 
olating this analysis to the MY model, where no conserved 
charge exists, he proposed that a physical definition of the 
number density could be given as Ns — [l/M)pBS, where 



Pbs represents those terms in the total energy density which 
have a Boltzmann-suppression factor. Moreover, he argued 
that non-Boltzmann-suppressed terms should be interpreted 
as corrections to the energy density of the x' s obtained after 
integrating out the heavy cj>'a at low temperature. Conclusive 
insight on the last point was achieved by Braaten and Jia 
[H], who worked in the effective lagrangian approach. They 
were able to show that the energy density of the x-fields actu- 
ally receives power-law corrections from the exchange of vir- 
tual <j>'s. However, the leading terms are only 0(X 2 T s /M 4 ). 
A 0(X 2 T 6 /M 2 ) term, which would correspond to the effect 
claimed by MY, was shown to be eliminated by a field redef- 
inition, thus unphysical. From the above discussion, we see 
that the problem of the definition of the number density for 
a heavy particle in a interacting theory is still an open issue. 
As discussed in |J the effective lagrangian approach cannot 
be of help in this respect, as it accounts for the effect of the 
heavy <f>'s on the light x' s , whereas what we need is a way to 
compute the effect of the x' s on the 0's, just the opposite. 

In this paper, we will show that it is possible to give a non- 
ambiguous definition of the number density of an interacting 
particle at any temperature and at any order in perturbation 
theory. In the case of the MY model of eq. (Ill), our definition 
exhibits the expected Boltzmann-suppression at low temper- 
ature. We will find it very convenient to work in the real-time 
formalism of finite temperature field theory (RTF) pi , which 
has the advantage of keeping T = and finite temperature 
effects well separated. Indeed, as we will discuss, the ambi- 
guities in the conventional definitions of the number density 
are mostly due to the T — O-renormalization procedure. Be- 
ing able to single out T = ultraviolet (UV) divergences 
turns then out to be a great advantage when one is interested 
in genuine thermal effects. Moreover, we will eventually be 
interested in using the defined number density in an out-of- 
equilibrium context, and in this respect the RTF, due to its 
close relationship with the Schwinger-Keldysh formalism, is 
clearly a much better starting point than the imaginary time 
formalism. 

The tree-level propagator for a real scalar particle in the 
RTF has the structure 

Do = D^ + Dj = zP[A ] + i(A - Ao)A/"(|go|)B , (6) 



where A = 



o = (l 2 ~ M ' 2 + (^)* £ )~ 1 is tne tree-level Feynman 
propagator, N the Bose-Einstein (Fermi-Dirac) distribution 
function, JV(x) = (exp(x/T) — (+)1) _1 and the two matrices 
P and B are defined as 



P[o(«)] = 



B 



a(q) (a - a*)(q)9(-q ) 

(a- a*)(q)G(q ) ~a*(q) 



1 1 
1 1 



(7) 

with 8(x) the Heaviside's step function. 

The full propagator has the same structure as the tree-level 
one, eq. (H), with A replaced by its all-order counterpart 



II (go, \q\) being the (renormalized) full self-energy at finite 
temperature. In the following, we will refer to iP[Ao] and 
iP[A] as the "T = 0" parts of the tree-level and full propa- 
gators, and to i(A - Aq)N (\q \)B and i(A - A*)J\f(\q \)B 
as the "thermal" parts. It is important to notice that what 
we call the "T = 0" part of the full propagator actually con- 
tains thermal corrections in the self-energy, so it is not the 
full T = propagator. 

Our discussion and results may be summarized as follows. 
We will start recalling how composite operator renormal- 
ization induces a mixing between the would-be number (or 
energy-density) operator for the 0-field and operators made 
up of x-fields (the "counterterms" in eq. (0)), which in general 
make a definition of a "pure" (jd-number quite cumbersome. 
These divergences, and the mixing, arise as the field opera- 
tors are evaluated at the same point. Diagrammatically, they 
arise when the ends of the two cj> legs are joined, provided both 
legs are T = propagators. It is inherently a T — effect, 
since if any of the two external legs is thermal, no divergence 
arises thanks to the UV cut-off provided by the Bose-Einstein 
function. Then, a thermal averaging procedure can be nat- 
urally defined, at least for quadratic operators, which is free 
from composite operator UV divergences. 

If we insist in defining the number density starting from the 
energy density p to t, we are led towards a different splitting of 
ptot with respect to that considered by MY. The p$, p x , and 
pint components into which we split ptot, are such that: i) each 
term is closed under operator mixing, so that the splitting is 
independent on the renormalization scale; it) no 0(A 2 T 6 /M 2 ) 
correction to any of the pieces in which ptot is split emerges; 
Hi) the </>-energy density is Boltzmann-suppressed. Finally, we 
will point out that a better definition for the number density, 
also exhibiting Boltzmann-suppression at T < M but valid 
in any temperature range, can be given by considering the 
thermal part of the full propagator in RTF. 

The paper is organized as follows. In sect. II the basic 
features of composite operator renormalization are illustrated 
for the case of the <fi 2 operator in the MY model. In sect. Ill 
a new definition of thermal average is given, which does not 
require composite operator renormalization. In sect. IV this 
definition is applied to the energy-momentum tensor in the 
MY model, showing that it leads to a splitting of the total 
energy density in <j>, x, and "interaction" pieces which differs 
to that considered by MY and is not plagued by the usual 
ambiguities induced by operator mixing. In sect V a defi- 
nition of the number density valid at any temperatures and 
extrapolable even to the massless limit is given. Moreover, 
it is shown how an extension of the MY model in which the 
heavy fields are complex and carry a conserved charge natu- 
rally leads to a Boltzmann-suppressed definition of the heavy 
particle number density. In sect. VI the case of a heavy par- 
ticle which is unstable already at T = is discussed, showing 
that power law corrections do indeed emerge in this case. Fi- 
nally, in sect. VII we discuss the relation of our work to that 
of refs. [0-p| and give our conclusions. 



A w (g) 
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II. RENORMALIZATION OF 4>l 



III. THERMAL AVERAGE WITHOUT MIXING 



Before considering the complete energy-momentum tensor, 
the analysis of the composite operator 0q (where the 0-label 
indicates the bare fields, masses and coupling constants) will 
be useful in illustrating the main points of renormalization 
and thermal averaging. We start by considering the T = 
case. Due to the local product of fields making up the op- 
erator, new UV divergences are induced when it is inserted 
in a Green function, which are generally not canceled by the 
lagrangian counterterms. The renormalization of these diver- 
gences then requires the introduction of new counterterms. 
Renormalized operators can be defined, which are generically 
expressible as linear combinations of all the bare operators 
of equal or lower canonical dimensionality (for a thorough 
discussion of composite operator renormalization, see for in- 
stance fhollllj] ) . In the model (DJ) the renormalization of 4>q 
then induces a mixing with Xo, an d the two renormalized op- 
erators, [4> 2 ] and [x 2 ] can be expressed as 



[x 2 






sz, 



4>X 

sz x 



(pt) 
xi 



(9) 



In minimal-subtraction scheme (MS) the bare parameters can 
be expressed in terms of the renormalized ones as 



Our ultimate goal is to build a definition of number density 
for <j> particles out of the bilinear operator <\> (or H^ as in the 
MY case). From the discussion in the previous section, we see 
that a major obstacle in this respect is the feature of operator 
mixing, since there seems to be an inevitable "pollution" by 
the x 2 operator (implied by SZ^ X ,SZ X ^, 7^ 0), preventing us 
to cleanly separate between the two particle species. More- 
over the /^-dependence of SZ^ X , SZ x< p makes any definition of 
number density renormalization-scale dependent. 

A diagrammatic analysis in the RTF can help us in iden- 
tifying the T = nature of such pollution, and to give a 
definition of thermal expectation value for (f> 2 which is free 
from it. 



Let us consider tr e 



-H/T 



[0 2 ]/tr« 



-H/T 



At tree-level in the 



RTF it is given by the two diagrams in fig. If 






(a) 



Aj,o — n Aj I 1 + 



y 4({A 3 };fl) 



1+ 
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m 



1 + 



sw_ 

M 2 
Sm 2 



= M' 1 + -^ 



P? V-fefdA,};/?) 



M 2 Z^, 

fc=i 



1 + 



E 



(D - 4) fc 
(D - 4) fc 1 



(10) 



where {Ai} = {A^, A x , A}, D is the space-time dimensionality, 
and the important point is that the a^'s, b^ 's and b™'s are 
mass, momentum, and temperature independent. Using eqs. 
(loft it is then possible to prove that, for any renormalized 
green function G the following relation holds 



M 



dG 

dM 2 



+ m 



dG 

dm 2 , 



M" 



dG 

dM 2 



+ m" 



dG 

dm 2 



showing that, being the right hand side manifestly finite, the 
insertion of the combination M 4> + rn Xo does not induce 
any UV divergence. The SZ^, 8Z XX renormalization con- 
stants can then be expressed in terms of the mass countert- 
erms as 



&Zsa — 



sz xx = 



SM 2 



M 2 
M 2 



SZ 



x4> 



SZ, 



4>x 



m A m A 

whereas the remaining counterterms SZ<f, x , SZ X 
puted as 

{4>l{x)x{v)x{z)) 



(11) 



can be com- 



8Z^ X = 



8Z x4 , = 



(x 2 (x)x(y)x(z)} 



{<ft(x)<Kv)4>(z)) 






+ (b) 



SZ^x 



SZ <)>x 





(c) 



FIG. 1. Diagrammatic expansion of ([0 ])t. 



where the continuos line represents the T = part of the (j> 
propagator, the barred continuos line the thermal part and 
the double line the sum of the two. The first diagram is 
quadratically divergent and requires a subtraction, whereas 
the second one, thanks to the Bose-Einstein function con- 
tained in D (q), is finite. At O(A) the diagrams of fig. fin con- 
tribute. There are also contributions O(X^) (such as higher 
order contributions containing also A x ) but their consider- 
ation is not necessary for our discussion, so we will limit 
ourselves to contributions containing only powers of A. The 
(quadratic) divergence due to the T = x-tadpole diagram 
(dashed line) is cancelled by the mass counterterm. The only 
remaining divergence comes from the first diagram, in which 



both the 4> lines, entering the cross are of T — type. It is 
the cancellation of this divergence which calls into play the 
mixing, via the renormalization constant SZfa appearing in 
the last diagram. At 0(A 2 ) (fig. nlc) the same happens, where 



We are now ready to define a number density for cf> which 
will not require any mixing with operators involving the x~ 
field and, as a welcome byproduct, will exhibit the expected 
Boltzmann-suppression at low temperature. 



is the self energy to 0(A 2 ) renormalized by lagrangian coun- 
terterms. The last diagram, containing a tree-level \ 2 opera- 
tor multiplied by a 0(X 2 ) renormalization constant is required 
in order to cancel the divergence coming from the integration 
over the momentum of <j) particles flowing to the cross with 
T — propagators (the first diagram in fig. lc). Note that, 



due to the general relation ^ . ._ £;j 



0, Eij being the 



self-energy with external thermal indices i and j, the con- 
tribution with both lines of thermal type also vanishes, once 
the summation over thermal indices is taken. The remaining 
loop divergences are cancelled by usual lagrangian countert- 
erms and therefore require no mixing. 

The only divergence which is not cancelled by neither la- 
grangian nor composite operators counterterms is the tree- 
level one. Then, one usually defines the thermal expectation 
value of <fi 2 as the expectation value of \d> 2 ] subtracted accord- 
ing to one of the two following ways |>lll2j , 



[<^]>T 



tr[0 2 



-H/T 



tre- H / T 



tr . 



W 2 ]|Q) 

<0|0> 

<O|[0 2 ]|O> 



tre-»/ T 



(0|0) 



(12) 



The definitions above are both finite, and thus a plausible 
choice. However, both of them include the contributions from 
the first and the last graphs in figs, fin and file, that is, they 
require mixing. 

The above analysis suggests a new definition of thermal 
average of (f> 2 , which represents the main point of this paper, 
namely 



2 \ 

/Th = 



A 



(2tt 

A 



(27T 



jD 1 (q) 

i(A-A*)(g)Af(M) 



(13) 



where D T (q) is any of the four component of the thermal part 
of the full propagator. 

Expanding A w as A w = A^ + A^IIA^ + ..., and 
recalling the relations between II(qo, \q\) and E(qo, \o\)ij P|, 

Ren(g) =ReEn(«), Im 11(g) = e(go) Im [En (g) + Ei 2 (g)] 

one can check that the perturbative expansion of fig. fij is re- 
produced, apart from the 'problematic' diagrams, i.e. the 
tree-level T = diagram, the higher order ones with both 
T = lines flowing into the cross, and the diagrams contain- 
ing composite operators counterterms, which are all omitted. 
The definition in (113) then automatically gets rid of all the 
unpleasant features of the conventional definitions in eq. (fl2[), 
namely, the need of an arbitrary subtraction and, more im- 
portantly, composite operator renormalization and mixing. 



IV. AVERAGING THE ENERGY-MOMENTUM 
TENSOR 



The renormalization of the energy-momentum tensor is dis- 
cussed in refs. [9, , 12] . Expressing it in terms of bare param- 
eters, the operator 

I*" = a"0o0"to + 9 M Xo9"xo - g^Co + ■ ■ ■ , 

is finite when inserted in a Green function, that is, it does 
not require extra counterterms besides those already present 
in the lagrangian. The dots in the formula above represent 
pole terms proportional to {d^d v - g^ 1 'd 2 )<j>l and to (<9 M <9" - 
g^ v d 2 )xx). Since we are interested in the thermal average of 
T^ v alone, translation invariance insures that such pole terms 
do not contribute. On the other hand, a subtraction of the 
T — contribution as in eqs. (113) is needed. Working with 
renormalized parameters, T^ v is expressed as 



rpfjt 



(14) 



where T££ n has the canonical form and T£% contains only the 
lagrangian counterterms. 

The finitness of T^ v means that the divergences induced 
by, e.g., the composite operator 4> 2 x 2 are canceled by the la- 
grangian counterterms and the particular combination of op- 
erators 4 , \ 4 y {d4>) 2 , ■ ■ •, appearing in T**". In general, if we 
split T Ml/ in parts, each part separately will require composite 
operator renormalization. Indeed, this is what is done by MY 
in ref. H. They split the total hamiltonian as 



rriU 



H t , 



Ha, + H x + _ffi n 



H c 



and compute thermal average of Hj, (given in eq. (2)). In the 
RTF it gives 



tr-H* 



-H/T 



lre- H ' T 



cTg 

(2tt)4 



[D J 



[2q 2 - (q 2 - M 2 )} 



+ D T ]i2(g)], 



(15) 

where, differently from eq. (|13|) , also the T = part of the 
full propagator appears. It is important to recall here what 
we have noticed after eq.(pl), namely that the T = full prop- 
agator contains thermal effects via the self-energy. Such con- 
tributions are not cancelled by purely T = subtractions 
like those defined in eq. (|12j). By perturbatively expanding 
[D T=0 ]i2 in (113) at 0(A 2 ) we find, apart from divergent con- 
tributions which require composite operator counterterms to 
be added to H^, the "famous" power-law contribution 



d 4 q 



[2% 



M 2 



(2 I )4 
§/(|^[2« 2 -(? 2 -M 2 )]0(-<zo) 



D T= °] 12(Q)! 



o(\ 2 ) 



(16) 



= div. + 



AJ Z ) Ren - (A 2 , + Ao 2 )e(q ) ton] 
1 A 2 T 6 
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where e(go) = 0(qo) — 9(—qo). The same contribution, with 
opposite sign, is found from the corresponding piece of the 
thermal average for H x . 

Paralleling our discussion in the previous section, we will 
instead define the energy density of the cf> field as 



rf 4 



(27^ 



[2g 2 -(g 2 -A/ 2 )][D T ] 12 (q). 



(17) 



and analogously for p' x . 
then be split as 

Ptot = (HtotjT = 



The complete energy density will 



trH t 



tote 



-H/T 



<0|#tot|0) 



tre-»/ T 

P4-+ Px + Pint , 



(0|0> 



(18) 



so that the contribution of eq. (nm enters p' int instead of p'^ . 
Notice that neither p'^ nor p' x are affected by the T = sub- 
traction. 

The definition in ( |l8|) exhibits three remarkable properties: 

i) The splitting in (|lq) is closed under operator mixing and 
composite operator renormalization. No composite operator 
counterterm is required to make p'^ or p' x separately finite and, 
since p to t is also finite, the same is true also for p' int . Thus, 
our definition of number density, differently from that consid- 
ered by MY, is independent on the renormalization scale p rl 
There is no "pollution" from p' x t o pL at any value of p. 

ii) None of the three pieces in ([18|) contains 0(X 2 T 6 /M 2 ) 
terms. Indeed, the contribution of eq. Qlq) and the analogous 
one for the x' s g° both into p' int and there, because of their 
opposite sign, cancel. 

Hi) p'^ is Boltzmann-suppressed for T <C M. Indeed, writing 



[D J 



fD J 



2(g) explicitly 



• («) = - 



2lmn(q)e{q ) 



(q 2 - M 2 -Re 11(g)) 2 + Im 11(g) 



r^(kol), 



(19) 



we see that, due to the Bose- Einstein function, non- 
Boltzmann-suppressed contributions to the integral in eq. ( |17J ) 
might only come from momenta \qo\ Ss T (<C M). But, for 
these values of momenta, it is the imaginary part of the self- 
energy which is Boltzmann-suppressed, as we show explicitly 
in the Appendix. 



V. A BETTER DEFINITION FOR N, 







Using the 0-energy density as given in eq. (|17J) , we can now 
define the <j> number density for T <C M in a way analogous 
to MY's definition, eq. (g), that is 

N * = M for T « M • 



However a better definition, valid at any value of T, and again 
based on the average of a quadratic operator, can be given as 



N<t 



d 4 q 

(2tt) 4 



W 



[D T ] 12(9). 



(20) 



As p'^, eq. (GOh exhibits the nice features of not requiring 
composite operator renormalization and of being Boltzmann- 
suppressed at low temperatures, but in addition it may be ex- 
tended to high temperatures or, equivalently, to the massless 
limit. The origin of eq. (M) can be understood considering 
the case of a complex scalar field, $, in case an exact (7(1) 
symmetry is imposed to the theory. In this case, a natural 
definition of number density involves the current operator 



Jv( x ) = g 



&(x) d„ $(x) - $(x) d,, &(x) 



In the RTF the thermal average of this operator corre- 
sponds to id Xfi D c {x - y)\ x=y where D c = [D]i 2 + [D] 2 i. 
The thermal part of the propagator for a complex scalar 
field differs from eq. (g) only in the replacement of Af(\qo\) 
with Af(q ) = 6(q )Af+(qo) + 0{-q o )Af-(q o ), where Af±(x) = 



(exp(/3(x =F M*! 



1) 



and p<f, is the $ chemical potential. 



Since the current is conserved, the evaluation at x — y does 
not require composite operator renormalization. Indeed, it 
can be checked explicitly that, after Fourier transforming, we 
get 



O'o(a0>T = 2 



go [D T ] 12(g). 



(2tt 



(21) 



where the T — contribution vanishes. Incidentally, this 
sheds new light on the conclusions of ref. M in which it was 
shown that in case a U(l) symmetry exists, a Boltzmann- 
suppressed number density can be defined. From what we 
saw above, it is so because in this case the current operator 
for $ does not mix with any other operator. In absence of 
mixing, terms like Qlq) are absent and the number density is 
automatically Boltzmann-suppressed. 

At tree-level we get the well-known result for the charge 
density 



= Uo(x))t = 



d 3 q 

j2lPp 



[A/-+H-A/1H] 



where w = (|gf + M 2 ) 1/2 . 

Now we can see how eq. (ECJ) can be derived in analogy to 
(Ejl). The number density for the real scalar field 4> should dif- 
fer from the previous quantity in two respects. First, particle 
and antiparticle contributions should be summed up, in order 
to account for their 'Majorana' nature. This is obtained by 
taking the modulus of |qo| inside (EOf). Second, a factor 1/2 is 
needed, in order to take into acccount the number of degrees 
of freedom. 



VI. THE CASE OF DECAYING HEAVY 
PARTICLES 



: Of course, at any finite order in perturbation theory we 
still have the usual p dependence induced by lagrangian 
counterterms. 



The discussion of the previous sections is straightforwardly 
generalizable to any model in which heavy particles annihi- 
late into lighter ones, both bosons and fermions. The ba- 
sic conclusions remain unalterated, the main reason being 



the Boltzmann-suppression of the imaginary part of the self- 
energy at small momenta (see the discussion at the end of 
Sect. IV and in the Appendix). 

The situation changes drastically if the heavy particle is 
unstable at T = 0. If, for definitness, it decays into two 
lighter particles of mass m < M/2 the propagator exhibits an 
imaginary part for 2m < po < M already at T = 0, so that 
there is no reason to expect that it is Boltzmann-suppressed 
for T < M. The heavy particle is actually a resonance, whose 
energy-momentum can take values much lower than the peak 
at p 2 — M 2 , of course paying the usual Breit-Wigner sup- 
pression far away from the peak. In this case, power-law 
contributions to the resonance number density do in general 
emerge, as we will briefly discuss. 

We will consider for definitness the annihilation model dis- 
cussed in ref. |7| and in the first of refs. ]3|, in which the 
interaction between the heavy $ and the light ^-bosons is 
due to the term 



Cin 



2 A 



(22) 



where we assume that p <C M. In this model, the imaginary 
part of the self-energy is oc fi 2 , so that, using ( |19| ) in (EOI) we 
get 



"-B^ 



(23) 



with To = h 2 /(32tvM) the T = decay rate on-shell. In ref. 
M a different behavior, iV$ ~ T 4 was claimed. Indeed, the 
following argument clarifies the physical origin of the power- 
law contribution and confirms the T 5 behavior. 

In the thermal bath, unstable $ particles are continu- 
ously produced in \~X annihilations. Being the $ a reso- 
nance, the production energy, u), can be much smaller than 
M. The produced $'s eventually decay with inverse lifetime 
T = fj, 2 J (32ttuj) - ToM/uj. Thus, the number density of $'s 
of energy around oj, n (uj), obeys the rate equation 



<in*(u;) 
~dt 



= 7-rn*(tj), 



where 7 is the rate per unit volume of the process 



XX 



$ — > everything 



(24) 



(25) 



The above equation leads to the equilibrium value 
n eq (uj) = - . 



N x av, where the inclusive 



The production rate is given by 7 
cross section a may be computed using the optical theorem as 
a ~ fi 2 ImA$/s, s being the square of the $ four-momentum, 
and A$ the full 4> propagator. 

Taking m <C y/s ~ T <C M, the initial \ particles are rela- 
tivists (i.e. N x ~ T 3 ) and we get a ~ p 4 /(M 4 T 2 ). Putting 
all together, we get a contribution to the total <f>-number den- 
sity from the region uo ~ T of order 






(26) 



The $'s with energies closer to the peak of the resonance 
live longer but their contribution to the number density is 



Boltzmann-suppressed, as two highly energetic x' s are re_ 
quired in the initial state to produce them. 

As the energy of the relevant $'s is typically 0(T), the 
energy density is 



p$ ~ O 



(Is. 

\ M 3 " 



T J + Boltzmann — suppressed , 



and pg>/p x vanishes as T 2 at low temperatures. 

Before concluding this section, we observe that the defini- 
tion of a sensible number density along the lines of eq. (pi), 
i.e. starting from the free hamiltonian, does not seem quite 
appropriate in this case. Indeed, as we have just discussed, 
the width of the $ plays a crucial role, allowing energy values 
much lower than the tree-level mass-shell. It seems that, if 
one wants to insist in looking for definitions like eq. (H) , a free 
Hamiltonian for quasi-particles, along the lines discussed for 
instance in ref. Hj, should rather be used. 



VII. CONCLUSIONS 

The definition of particle number density in an interacting 
theory is a delicate matter. The necessity of giving a mean- 
ing to divergent composite operators calls into play operator 
mixing, so that a separation between different particle species 
turns out, in general, to be a renormalization scale dependent 
procedure. 

Our guideline in this paper was precisely to find a definition 
of number density free from this problem. We discussed how 
composite operator renormalization is inherently a zero tem- 
perature effect which can be isolated from the physically rel- 
evant finite temperature effects. We showed that in the RTF 
a proper definition of the energy and number densities comes 
out quite naturally and can be systematically computed at 
any order in perturbation theory. 

In the framework of the MY annihilation model, our defini- 
tion of thermal average for bilinear operators leads to a differ- 
ent splitting of the total energy density with respect to that 
considered by MY. Each piece is finite, i.e. it does not require 
extra counterterms with respect to those already present in 
the lagrangian. Moreover, as a welcome byproduct, we obtain 
a Boltzmann-suppressed energy density for the heavy particle 
<j>. Finally, no power-suppressed terms of the type found by 
MY appear. In this respect, we agree with refs. W-M]; the cru- 
cial point is how the interaction is treated. In our discussion, 
we concluded that diverging contributions like eq. (hd) should 
naturally go into the 'interaction' part of the total energy, to 
make it finite. If this is done consistently, 0(\ 2 T 6 /M 2 ) terms 
cancel out in p' int . This is consistent with the conclusion by 
Braaten and Jia M that these terms are unphysical. 

Our definition of the (^-energy density is close in spirit to 
that proposed by Srednicki in H. However the detail is dif- 
ferent. Namely, in Srednicki's definition, all the Boltzmann- 
suppressed terms should be included into the 4> energy density, 



the contribution depicted in fig'M 




FIG. 2. A diagram giving a Boltzmann-suppressed contri- 
bution to p[ nt 

The interesting question, now, is how to obtain quantum ki- 
netic equations giving our number densities in the equilibrium 
limit. From the present investigation it appears that the RTF 
may provide useful physical insights in this regard. 
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where C = j^^g e(go), wi = |g - p - fe|, ^2 = |p1 and 

CJ3 = v|fc| 2 + M 2 . Since Imll is an odd function of go we 
can restrict to go > 0. Imll can then be decomposed as: 

Imn = IhiiId + ImIL4 + lraYl LD 
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APPENDIX 

The imaginary part of the <f> self-energy in the MY model 
emerges at two-loops by the diagram in fig. til 




FIG. 3. Two-loop contribution to the self energy . 

From eqs. (fnl), and (h9|) we know that the contribution to 
the integral defining p<L from energies |go| > T is Boltzmann- 
suppressed. Non Boltzmann-suppressed contributions might 
only come from |go| 5s T(<C M), but we will show at the end of 
this appendix that this is not the case as Im II is Boltzmann- 
suppressed in this regime. 

According to the standard rules of the RTF M we have for 
the imaginary part: 



ImII(g)= -i e (g )( e ^°-l)Ei 2 (g), 



where 



A 2 

Si 2 (g) = — - 



^j^D 12 (p)D 12 ( q -k-p)G 12 (k), 



and 



iD 12 {p) = 2ir($(-po) +Af{\po\))5(p 2 ) 
iGrz{k) = 27r((9(-fco) + N{\k \))5{k 2 



M" 



(27) 



the x an d 4> particle propagators respectively. Integrating 
with respect to fco and po we obtain: 



Imn(g) = C 



(fpcPk 



{ 



LU1U2U3 

[5(g — (wi +LU2+ W3)) — 3(qo + +ui + W 2 + 0J3) 



where: 



Imlli 



C 



d 3 pd 3 k 



5 (go — (wi +lo 2 + con)) 



[(l +M)(i + A£)(i +A/ 3 ) -MA/iA/i] 



ImU A = C 



d 3 pd 3 k 



<5(go — wi - L02 + L03) 



WlW 2 U>3 

[( f + M ) ( 1 + A/2 )Af 3 - M A/i ( 1 + A/3 )] 



ImlL 



<5(go —CJ1+CJ2 —L03) 



2c (**a 

[(1 + M)Af 2 (l +A4) - M(l + A/iJA/k] + 



-c 



d 3 ptf 3 fc 



(5(go + wi + c;2 — W3) 



[(l+A/"i)(l+A/2)A/3-MAr 2 (H-A/3)] 



(29) 



with A/i = N(uJi). IhiIId corresponds to the decay of an off- 
shell 4> into two on-shell x' s and one on-shell <f> (minus the 
inverse process) (fig. kla), IuiITa to the annihilation process, 
where an off-shell </> annihilates with an on-shell <f> to produce 
two x's (fig. Mb), and Im TIld corresponds to Landau damping 
of the 4> induced by the scattering with a on-shell \ (fig. Wc). 



(a) 




(b) 



(c) 



FIG. 4. a) Decay, b) Annihilation, c) Landau damping . 
The computation proceeds as in |13j and, for q — 0, we get 

lmn D = C"0(g o - M) I dvf(v)f(uJ - v) ■ 



fiu^fiu' + v-u)' 

+ 00 



hi 
Imni = C"[% - M) 

[(/(«) + l)/(<J+r)l,, 

IhiIIlzj = 



[7] S. Matsumoto and M. Yoshim ura, Phys. Rev. D61 



dv + 9{M - q ) / dv] 

■J a 



C'e(M-qo) I dvf(v)(l + f(v-u))- 



(l + f(co 6 ))f(^ 5 +v-u;) 

^ + 00 



[8] 
[9] 

[10] 

[11] 
[12] 
[13] 



123508 (2000) Jhcp-ph/9910393 



H.A. Weldon, |hep-ph/9809330 , 

N. P. Landsman and C. G. van Weert, Phys. Rept. 145, 

141 (1987). 

J. C. Collins, Renormalization , (Cambridge University 

Press, 1984), Chap. 6. 

L. S. Bro wn, Annals Phy s. 126, 135 (1980). 

P. Jizba, 



hep-th/9801197 



E. Wang and U. Heinz, Phys. Rev. D53, 899 (1996) [hep- 
ph/ 9509333]. 



2C"[%o-M) 
[f{v)f{u - v) In 



dv + 6(M-q ) 

(1 + /(^ 6 )) ■ 

(1 +/ („_ W+W 6))J 



dv] ■ 
(30) 



with 



C" 



X 2 



128tt : 



;e(«o)(e 



/3<?o 



1)£. 



and 



M 

go 
T' 



ql + M 2 



2g T 



0^3 



1,2 1 90 

W = J — 



W3T* 3,4 

- (J = 



r 



1 ^3 
2 



go T& 



T 



! go + u 3 T j» 

2 rp 



(31) 



From the formulae above it is easy to check that the ImiT is 
Boltzmann-suppressed for jgo| Si T(<C M). Physically, it can 
be understood by looking at fig.^. For jgoj Ss T(<C M) only 
the annihilation and the Landau damping contribute. In the 
former case (fig. Wo) the on-shell <j> particle in the initial state 
comes from the heat-bath and then carries a N factor. In the 
Landau damping case (fig. fie), the energy to create the on- 
shell <j} has to be provided by the x from the heat bath, so 
that a Boltzmann-suppressed N has to be payed in this case 
too. 
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